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Abstract 

The dependence of the transmission coefficient on the deformation, the collective 
rotation and excitation energy of the compound nucleus emitting light particles is 
introduced in the framework of Weifikopf's evaporation theory. The competition 
between fission and particle evaporation is treated by a Langevin equation for the 
fission variable coupled to the emission process. Detailed calculations are presented 
on the decay of different Gd and Yb isotopes at an excitation energy of about 
250 MeV. These calculations demonstrate the importance of the effects of nuclear 
deformation and of the initial spin distribution on the evaporation. 
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1 Introduction 



The properties of excited nuclei in thermal equilibrium are of great physical interest. 
We are interested in excitation energies below some (400-500) MeV, where we may still 
describe the nucleus as a system of neutrons and protons which interact by effective forces. 
All the information on the physical state of hot nuclei is then to be obtained from a careful 
study of their decay by emission of neutrons, protons, and a-particles and their decay by 
fission. The emission of photons will be neglected, as we consider energies far above the 
thresholds for particle emission. 

The basis of our approach [0 is a description of the particle emission as a purely 
statistical process, as given by the Weifikopf theory M , and the nuclear fission as a trans- 
port process 0. The importance of the non-statistical aspects of the fission process in 
this context was recognized by Grange and Weidenmiiller ||. Our approach and also the 
ones of Frobrich, Abe and Carjan [[5], ^|, [7j are based on the premises of their work. 

Our efforts are dedicated to a stepwise improvement of this general theory in order 
to reach a quantitatively reliable description of the decay of the compound nuclei. 

In the present paper, we give a more careful study of how the evaporation proba- 
bilities depend on the deformation of the nucleus, on its excitation and on its collection 
rotation. In the standard form of the evaporation theory @, the deformation of the de- 
caying nucleus enters only through the level densities of the initial and final nucleus. We 
take into account that the transmission coefficient also depends on the deformation and, 
to a smaller degree, on the collective rotation of the nucleus. This will be explained in 
Chapter 2. 

Our research is encouraged by the increasing amount of careful experimental studies 
of the evaporation of light particles from excited nuclei and of the concomitant decay by 
fission ||. The results of our calculations are presented and compared with experimental 
work, especially of that of Ref. [[| in Chapter 3. 

Finally, we summarize in Chapter 4 our findings and point out the direction of our 
future work. 



2 Description of Fission Dynamics Including 

Evaporation 

The fission process of hot and rotating nuclei should be described within a statis- 
tical model which takes into account the effect of energy dissipation due to the coupling 
to the internal degrees of freedom modelled by a friction force which generates diffusion. 
The corresponding transport equation of the Fokker-Planck (FPE) type was originally 
proposed for nuclear fission by Kramers in Ref. O and later on adapted for heavy ions 



physics in Ref. [ 10 1 - [12]. It is not easy to solve this equation exactly when the num- 



ber of collective coordinates and conjugate momenta is larger than 2. An effective but 
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approximate method of solving the FPE based on the moment expansion was used in 
pnH . Unfortunately the probability distributions in the multidimensional space of collec- 



tive coordinates is usually far from a gaussian form and the moment expansion method 



fails | 13fl . It is important to recall here that the set of differential equations which one 
obtains in the moment expansion method for the average coordinates and momenta is 
simply the set of classical equations of motion with friction as e.g. those used by Blocki 
and co-workers [J14jj. A more efficient way of solving the transport equations based on 



Monte Carlo method has been proposed in Ref. JK]. It was also shown in |I| that the 
Monte Carlo calculation combined with the local moment expansion method transforms 
the FPE into a set of equations of motion containing random forces which is equivalent 
to a set of coupled Langevin equations. The Langevin equation has already been used in 
TT| in order to describe heavy-ion collisions and fission |HJ . This idea is very fruitful and 



allows us to describe many phenomena occurring in heavy-ion physics and fission due to 
statistical fluctuations (see e.g. the review article 

In Chapter 2 we will present the physical basis of our model and explain how we 
couple the emission of light particles to the dynamical evolution of the system from the 
compound nucleus state to the scission point. This dynamics is governed in our descrip- 
tion by the Langevin equation which we will present in Sect. 2.1, before explaining in 
Sect. 2.2 how particle emission is incorporated into this description. In Sect. 2.3 and 2.4 
we deal with the initial conditions and the integration of the Langevin equation. 

2.1 Equations of Motion 

We need to describe the time evolution of the nuclear system from an initial state which 
corresponds to a compound-nucleus at high excitation and angular velocity but usually 
created close to spherical symmetry to a final state in which a large amount of excitation 
and angular momentum will have been dissipated. For such a description we use a single 
collective coordinate q = p C m/Ro which measures the distance between the centers of mass 
of the two halves of the fissioning nucleus in units of the radius Ro of the corresponding 
spherical nucleus. For the time being we restrict our description to symmetric fission. 
The description of asymmetric fission would necessitate the introduction of an additional 
collective coordinate describing the mass asymmetry of the two emerging fission fragments. 
If p designates the conjugate momentum associated with the collective coordinate q we 
obtain the following equations of motion describing the time evolution of the fissioning 
nucleus: 

d q _ p m 
dt M{q) ' 1 ' 

dp _ 1/ p \*dM dV 7 (g) F( ) (2) 

dt 2\M{q)J dq dq M{q) V lU ' U 

In these equations M(q) represents the collective mass, V(q) the potential and 7(g) the 
friction coefficient [^l]]. The potential V is calculated as the difference between the 
Helmholtz free energies of the deformed and spherical nucleus. The friction coefficient 
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7 is calculated in the framework of the wall and window friction model |TB, [19]]. The 
explicit expressions of these quantities, as obtained in the framework of the Trentalange- 
Koonin-Sierk shape parametrization [2D|, have been presented and discussed in detail in 
Ref. ||21|| . The quantity F^(t) designates the random Langevin force which couples the 
collective dynamics to the intrinsic degrees of freedom. This Langevin force is chosen 
to be a gaussian random variable with zero mean value. In practice eqs. ([!]) - (0) are 
discretized by introducing a sufficiently small time step r. Then the equations take the 
form 

P 



q(t + r)-q(t) 



M(q) 



p(t + r)-p(t)^-{ 



V 



d M 



d V 



7(g) 

M{q)' dq dq ' M(q) 



r — 



t — 



pr + JD(q)f L (i 



(3) 



(4) 



Here — V^r] where r] is a gaussian distributed random number with <rj>= and 

< rj 2 >= 2 and where the brackets represent ensemble averages. The diffusion coefficient 
D(q) appearing in eq. @ is related to the friction coefficient 7(g) through the Einstein 
relation 

D(q) = 7 (g) T . (5) 

This relation holds in linear response theory as a high-temperature approximation of the 
dissipation-fluctuation theorem. In the applications presented below the nuclear temper- 
ature will always be sufficiently high so that eq. (BJ) is approximately verified. To speak at 
all about a nuclear temperature already supposes that the excitation of the nuclear sys- 
tem is shared to equal parts by all of the nucleonic degrees of freedom, so that statistical 
models apply. To simplify the description we consider the excited nucleus in a statistical 
description as member of a grand-canonical ensemble which can be characterized by 
a temperature T. For such a description to make sense, we need to assume that the time 
scale which governs the fission dynamics is sufficiently long compared to the one which 
determines the internal equilibration of the nuclear excitation among all the nucleons. If 
this is the case, we can assume that the system can be considered as being continuously 
at equilibrium. In this framework the temperature T is simply a measure of the nuclear 
excitation energy E* and related to the the latter through the usual Fermi gas relation 
E* = a(q) T 2 where a(q) is the level density parameter of the considered nucleus at a 
nuclear deformation characterized by q. The excitation energy itself is determined by the 
conservation of the total energy as will be discussed in section 2.4 below. 



2.2 Particle Emission 

The process of light particle emission from a compound nucleus is governed by the emission 
rate at which a particle of type v (neutrons, protons and a particles are considered 
here) is emitted at an energy in the range [e a — 4^ , e a + before the compound 
nucleus eventually undergoes fission. 
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Several theoretical approaches have been proposed in order to describe the emission 
from a deformed, highly excited and rotating nucleus Jl|, || [7|]. The method which we 
present below is close in spirit to the one used in |TJ , but in our description the widths 
for particle emission will in addition depend on the deformation and angular momentum 
of the compound nucleus p2[ . 

According to Weisskopf 's conventional evaporation theory || the partial decay rate 
r%P(E* , L) for emission of a light particle of type v with energy e a and orbital angular 
momentum E@ from a compound nucleus with excitation energy E* rotating with an 
angular momentum L can be written as 



L+lf, e a +Ae a /2 



where 



T a /{E\L)= ^ J w v (e,£ /3 ; X )p R (E* R ,L R )de , (6) 

P\ ' I L R =\L-e \ g a — Ae Q / 2 

p(E*, L) = (2L + 1)(— f'^a (7) 

is the level density in the emitting nucleus. The level density p R in the residual (daughter) 
nucleus with the excitation energy E R and the angular momentum L R is obtained in the 
same way. Both these quantities also depend on the mass and charge number and on 
the nuclear deformation. The quantities J and a represent the moment of inertia and 
the level density parameter of the compound nucleus and S v the intrinsic spin of the 
emitted particle. w u (e,£; x) is the transmission coefficient for emitting a particle of type 
v, with energy e and angular momentum I from the deformed compound nucleus. The 
parameter \ in the argument list of w v stands for all quantities not explicitly mentioned 
here, such as the mass and charge number, the nuclear deformation, the direction in 
space in which the particle v is emitted. Proceeding in this way would, however, leads to 
hardly tractable numerical problems, since in the Langevin formalism which we endeavour 
here, we need to follow the dynamics of the fissioning nucleus plus evaporation for a very 
large number of trajectories (of the order of 10 6 ). We therefore use a simplified procedure 
which introduces a transmission coefficient w v (e, £; x) obtained by a double averaging over 
the different emission directions and over the whole surface of the deformed compound 
nucleus. A detailed description of how w u (e,£;x) is calculated in the framework of the 
Hill- Wheeler approximation [^3J and how the averaging procedure is carried out is given 
in Appendix A. With this averaged transmission coefficient w the width for emission 
of a particle of type v and energy e a reads: 



e a +Ae a /2 



•I V _L f 

K(E *> L) = 2nhp(E*) J <"^x)PR(E R )de , (8) 

e a -Ae a /2 

where 

HK ' 2L + 1 
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is the angular momentum independent part of the density (|?P and a similar relation 
holds for pr. The effective transmission coefficient in (§) is obtained by performing a 
summation over all allowed angular momenta of the emitted particle and those of the 
daughter nucleus: 

1 imax L +*9 

< ff (e; X) = 7Tj—7 E E (2£* + l)^(e, lp\ x) , (9) 

Zlj + 1 tp=0 L R =\L-£fs\ 

here £ ma x is the maximal angular momentum available for the particle having the energy 
e. As already mentioned the emission width also depends on the mass and charge 
number A and Z as well as on the deformation of the compound nucleus. 

Once the emission widths known, one can establish the emission algorithm which 
decides at each time step [t, t+r] along each of the trajectories if a particle is emitted from 
the compound nucleus. It is the value of the emission width which ultimately have to 
decide which of the light particles is emitted and at which energy. Since represents 
the rate at which a particle of type v and energy e a is emitted, the total emission rate for 
a particle of given type, irrespectively of the energy at which it is emitted, is given by 

n q _|_ -1 em r ax n 

UE*,L) = J wt ff (e; X )p R (E R )de = ^ V (E*,L) , (10) 

where we have replaced the upper integration limit by some large enough constant e max 
(which will, a priori, be different for the different particles) and at which the probabilities 
for particle emission will essentially have vanished, n is the number of energy bins of width 
Ae a in the interval [0, e max }. The emission rate T for emission of any kind of particles is 
the sum of the T u : 

r = r n + r p + r Q . (11) 

First of all we need to decide whether a particle is emitted at all in the given time interval 
[t, t + r]. The probability for emitting any particle is given, for a small enough time step 
r, by 

P(r) = l-e~ TT ^Tt . (12) 

One then draws a random number rji in the interval [0, 1]. If r]i<P(r) a light particle 
is emitted. If the time step r is chosen sufficiently small, the probability for emitting 
a particle will be small. In this way we guarantee that in a time interval at most one 
particle is emitted and we avoid to consider the emission of more than one particle in 
each time interval. 

In the case that a particle is emitted one needs to decide next of which type this 
particle is. To this purpose one draws a second random number r/ 2 in the interval [0, 1] 
and determines the localisation of 1] 2 with respect to the covering of this interval by the 
three bins T n /r, T p /T, and T a /T. Depending on the bin in which 7? 2 is located, a neutron, 
or a proton or an a particle is emitted. 
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We still have to determine the energy with which this particle is emitted, and we 
accomplish this in the following way. We introduce the quantity 

which represents the probability that a particle of type v is emitted with an energy smaller 
than e a . The quantity covers the interval [0, 1]. Subdividing the interval [0, 1] in a certain 
number of equal bins, one decides with which energy the light particle is emitted by 
drawing a third random number 7/3 in the interval [0, 1]. Inverting the function in eq. (|T3|) 
IT" 1 (773) will fall in one energy bin 

K\vs)e[e a -^,e a + ^} (14) 

of the total interval [0, e max ] and thus decide on the energy with which the particle is 
emitted. 



2.3 Initial Conditions 

The description of the fission process including particle emission starts, in principle, from 
an initial state corresponding to a spherical or deformed compound nucleus whose shape 
is characterized by the collective coordinate go, the corresponding conjugate initial mo- 



mentum p , the excitation energy Eq (with corresponding temperature T = J _Eo/ a (<?o)) 
and the initial total angular momentum L . In practice the excitation energy E* may 
be known approximately from the experiment. This is, however, not the case for g , p 
and the angular momentum L . In the calculations presented and discussed in Sect. 4 we 
choose go corresponding to a spherical nucleus and draw the conjugate momentum p for 
each trajectory from a normalized gaussian distribution 

PM = 7«m eM -A ] - (15) 

Here m is the collective inertia at the deformation go. It is in principle experimentally 
possible to determine the distribution of initial angular momenta for instance by measur- 
ing the 7 multiplicities |24]]. As these distributions are rarely available, we chose L to 



have a fixed value but we shall discuss below the dependence of the final results on the 
choice of L Q . 



2.4 Integration of the Langevin Equation 

Once the initial conditions are fixed one can integrate the system of equations of motion 
(1-2) using their finite difference version (3-4). At each time step [t,t + r] one draws a 
random number 77 from a gaussian distribution which defines the fluctuating force in eq. 
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(4) and thus generates a trajectory in the variable q. Simultaneously one asks, as just 
explained, whether a light particle is emitted and decides, eventually, of its type v and 
energy e a . 

At the beginning of the whole process the system has a fixed total energy E to t which 
is the sum of the collective kinetic and potential energy 

2 

EcoU = + V(q) 

and a collective rotation energy 

E _L> 

^rot 



2 J 



corresponding to a rigid deformed rotator of moment of inertia J |23, [26| . The excitation 
energy E* at time (t + r) and thereby the temperature of the nucleus is redetermined at 
each time step through the conservation of the total energy 

Etot = E^a + Ef ot + E* = E { / 1 + E<£ +E} + e a + B u + E recoU , (16) 

where the indices i or f refer to the initial (time t) or final state (time t + r). B u is the 
binding energy of the emitted particle (which is zero unless an a particle is emitted), e a 
is its kinetic energy after emission and E recoa is the recoil energy of the nucleus due to the 
emission process. These three quantities enter the equation only when a light particle is 
emitted in the time interval [t, t + t}. Since the recoil energy is in any case very small, it 
will be neglected. Additionally it is assumed that the deformation as well as the angular 
frequency of the nucleus is not changed due to the particle emission. 

It is evident that each emission of a light particle carries away excitation energy and 
angular momentum and thereby increases the height of the fission barrier of the residual 
nucleus which, in turn, renders the fission event less and less probable. 

For each choice of the initial conditions one generates in this way a separate trajec- 
tory which is followed through in time. Emitted particles and their energies are registered. 
Each trajectory can either lead to fission if it overcomes the fission barrier and continues 
on to the scission point, or can end up as a rather cold compound nucleus if too much 
energy has been lost to make the crossing of the fission barrier possible. Sampling the 
total number of trajectories Nfi SS which have led to fission defines the fission cross section 
as <7fiss = Nfi SS /N with N equal to the total number of trajectories. 

Let us finally mention that it is important to take into account the fact that the a is a 
composite particle and that its emission therefore presupposes its existence in the nucleus 
prior to emission. This is in principle quantified by the introduction of a preformation 
factor whose value can in principle be determined ||27|| . Since this quantity is, however, 



very poorly known, we fix it in the following by multiplying the emission width r a by a 
factor f a (0 < f a < 1). The choice of will be discussed in the following. 
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3 Results 



We study the decay of compound nuclei of several isotopes of 64 Gd and of 70 Yb at various 
excitation energies ranging from 150 MeV to 300 MeV. We selected these nuclei because 
a careful experimental investigation of their decays is available [7] and because they are 
known to exhibit large ground state deformations. Consequently, it is of special interest 
to investigate the influence of deformation on the emission of n, p and a particles from 
excited states of these nuclei. 

The results concern the following physical aspects: 

i The dependence of the emission probabilities of neutron, protons, and a-particles on 
the deformation, the excitation energy and the angular momentum of the emitting 
nucleus (Figs. 1 to 9). 

ii The dependence of the emission probabilities n^(e), eq. (]13|), of a particle of given 
type v on the kinetic energy e (Figs. 10 to 15). 

iii The dependence, at different temperatures T, of the number of fission events and 
light-particle multiplicities on the time t (Figs. 16 to 18). Although these functions 
are not measurable, they are important for the physical understanding of the decay, 
especially the transient time phenomenon [4] and its dependence on the temperature. 

iv The spectral distribution of the emitted particles (Figs. 19 to 21). 

v The dependence of the multiplicities, fission cross section and barriers on the initial 
angular momentum of the nucleus (Figs. 22 to 25). 

vi The influence of the friction forces on the multiplicities of the emitted particles (Fig. 
26). 

We will now give a detailed description of these points. 

3.1 Dependence of the emission probabilities on the nuclear 
deformation 

In Figs. 1-3, the emission width T u , eq. (0)) f° r n , P, an d a-particles is shown as 
a function of the deformation p cm /Ro of the emitting source. The initial ensemble of 
decaying nuclei consists of 7 Q°Yb 90 nuclei at an excitation energies of 50 MeV, 150 MeV 
and 250 MeV respectively, and with a rotational angular momentum L = 40h. All the 
emission rates are seen to grow as a function of increasing deformation. This trend can be 
easily understood as for increasing deformation the transmission occurs through a larger 
surface. This effect has already been observed for all three particles PSJ. We notice, 
however, that the emission width for a-particles increases more steeply than the one for 
n and p for the two lower excitation energies E* = 50 MeV (Fig. 1) and E* = 150 
MeV (Fig. 2). This is due to the fact that, as the nucleus is elongated, the barrier 
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height for charged particles is reduced in the section of the surface which is farther away 
from the nuclear center and increased in the section which are closer to the center. The 
section where the barrier is diminished is larger than the section where it is increased. 
Consequently, the emission rate for charged particles increases faster than the one for 
neutrons. 

In Figs. 4, 5, and 6 the deformation dependences of the emission rates for n, p, and 
a particles are shown for different values of the rotational angular momentum, varying 
between L = Oh and 60 h. For n and p the influence of the rotation on the emission rate 
is seen to be negligible, whereas the emission rate of a particles grows by 10%-20% as the 
rotational angular momentum increases from h to 60 h. It is clear that the centrifugal 
force, which helps to overcome the barrier, is largest for the a-particle. Furthermore, the 
centrifugal effects matter the more the larger the barrier. The barrier for the a-particle 
is larger than the one for n and p. 

Let us point out that the rotational angular momentum of the nucleus has a very 
noticeable influence on the height of the fission barrier which decreases as a function of 
increasing angular momentum. Thus, at high angular momentum, nuclear fission can 
compete more effectively with evaporation. 

In Figs. 7-9, the deformation dependent emission width for n, p and a is shown 
for two different isotopes of both 64Gd and 70 Yb. The results can be easily understood: 
for given proton number, the emission width T n for neutrons is the larger the larger the 
neutron number (see Fig. 7). On the other hand, for given neutron number, the emission 
width T p for protons is the larger the larger the proton number (see Fig. 8). At given 
proton number, the emission width T a , for a-particles decreases with increasing neutron 
surplus whereas, at given neutron number, r a grows with increasing proton number (see 
Fig. 9). This is a simple consequence of the fact that the a-particle contains an equal 
number of neutrons and protons. 



3.2 Behaviour of the probability II^(e) 

An interesting quantity is the probability ILy(e) to emit a particle of type v with an energy 



smaller than e. Its definition in terms of the emission rate T u (e) is given in eq. ([HJ). The 
emission width Y v depends also on the shape of the emitting nucleus, and so does the 
integrated probability Il^e). 

In Figs. 10 to 12, the probability Tl v (e) is shown for the emission of neutrons, 
protons, and a-particles. In each of the figures, n„(e) is plotted separately for the case 
that the emitting nucleus has a spherical shape, the shape corresponding to the saddle and 
the scission point. In all cases the fissioning nucleus is 160 Yb at an initial excitation energy 
E* = 250 MeV and an initial angular momentum L = 40h. All curves show a monotonous 
increase of U u as a function of the energy e starting from a minimal energy e which is 
zero for neutrons, but finite for protons and ct-particles due to the acceleration of the 
charged particles in the Coulomb field of the residual nucleus. For neutron emission, the 
dependence of II n (e) on the deformation of the emitting nucleus is very small, whereas for 
protons and a-particles the curves ILy(e) are shifted to somewhat higher energies for the 
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spherical source. The physical interpretation is that the Coulomb barrier is larger for the 
spherical nucleus than for the deformed nucleus and consequently the threshold for proton 
or a-emission is at a higher energy than for neutrons. We note that the scission shape 
for a nucleus like 160 Yb consists essentially of two tangent, slightly deformed fragments of 
half the charge. 

In Figs. 13 to 15 the probability Il^e) is shown for neutrons, protons, and a-particle 
emission for 3 different excitation energies of the initial compound nucleus 7o°Yb. The 
shape of the emitting nucleus is chosen here as corresponding to the saddle point and its 
initial angular momentum is fixed at L = 40 h. It is seen that the rise of the probability 
Ily as a function of e is the steeper the smaller the excitation energy. Such a behaviour 
can be easily understood since the range of energies of the emitted particles must rise with 
excitation energy. We have also investigated the dependence of the functions n^e) on 
the neutron and proton numbers of the initial ensemble of compound nuclei, comparing 
the emission from the 4 nuclei g4 4 Gd, g| 4 Gd, 7 o°Yb, and 7 o°Yb, in each case for the same 
excitation energy (150 MeV), the same angular momentum (40 h) and the same shape 
(saddle point). The results are indistinguishable for the emission of neutrons and only 
very slightly dependent on the nature of the emitting nucleus for the emission of protons 
and a-particles. 



3.3 Time dependence of the multiplicities and the number of 
fission events 

The dependence of the number of decays of a given type on the time which elapses 
starting with the formation of the compound nucleus is unfortunately not measurable. 
Nevertheless, we think that it is interesting to exhibit this dependence for a few cases 
since it enables us to gain insight into the dynamical mechanism. 

In Fig. 16 we show the number of fission events as a function of the time t on a 
logarithmic scale. This result was obtained with the light-particle evaporation channels 
turned off. The initial compound nucleus is ^Yb with an initial angular momentum 
L = 40 h. The 3 curves in Fig. 16 correspond to 3 different initial temperatures resulting 
in 3 different initial fission barrier heights Ub- It is seen that the transient time increases 
with decreasing excitation energy, as one expects. It should be noted that the functions 
Nss S (t) obtained for T = 4 MeV look quite similar to those for T = 5 MeV and are just 
shifted along the log(t) axis. Furthermore, if to is the time where half of the final number 
iVfi ss (oo) of fission events have occurred 

AWto) = ^ fiss (oo) (17) 
we may approximate iVfi ss for times close to to by a linear function of log(t) 

iVfissOO « N &ss (t ) + K log- (18) 

to 

As one infers from Fig. 16, the dependence on the initial temperature is mainly contained 
in the quantity to whereas the factor k,q in (ITS) is almost the same for T = 5 MeV and 
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T = 4 MeV. The approximate time dependence ( |i8| ) which is valid during part of the 
transient time interval is seen to be totally different from the one in the Kramers regime 
which is valid in the cases when the fission barrier is much higher than the temperature 
of the fissioning nucleus. 

In Fig. 17 we show the multiplicity (number of emitted particles in coincidence 
with fission and per compound nucleus) for the emission of neutrons, protons and a- 
particles as a function of time for the decay of a compound nucleus 7o°Yb at an initial 
excitation energy of 293 MeV and an initial angular momentum of 45 h. It is seen that the 
emission of p and a-particles ceases already after some 3-1CT 21 sec whereas the number of 
emitted neutrons still increases. Again this is easily understood as a result of the different 
thresholds for charged and uncharged particles. 

In Fig. 18 the fraction Nr ss /N of nuclei undergoing fission is shown for the initial 
compound nucleus 7 o°Yb as a function of time. Now, contrary to the results presented 
in Fig. 16, the emission of light particles is taken into account. The initial excitation 
energy of the nucleus is E* = 293 MeV and the initial angular momentum L = 45 h is 
assumed. The well-known Kramers approximation holds whenever there is an approx- 
imately constant current across the fission barrier which implies a linear dependence of 
the number of fission events upon time. It is seen from the figure that there is in fact 
no clearly distinguished section with a linear time-dependence. At best one could replace 
the function by a straight line in the interval between 8 and 12-10 -21 sec. The "transient 
time" , i.e. the time needed for the build-up of an approximately constant fission current 
is seen to be ~ 8-10 -21 sec. This is the main part of the time available for fission events. 
It is thus clear that the use of the Kramers approximation or of the still simpler statistical 
transition state hypothesis for the whole time interval would lead to wrong results. 

Furthermore the very concept of a constant current across the fission barrier, which 
is the pre-requisite of the Kramers approximation could loose its validity because of the 
existence of light particle emission which implies that the excitation energy and the an- 
gular momentum changes for those compound nuclei which emit a light particle prior to 
fission. 

3.4 Spectral distribution of the emitted particles 

In Fig. 19, we present the probability per energy unit of neutron, proton, and a-particle 
emission as a function of the energy of the emitted particle. This figure is made on the 
basis of 10 6 trajectories. The solid curves show the spectral distribution in coincidence 
with fission, i.e. for particles emitted from nuclei which subsequently undergo fission. 
The dashed curves represent the spectral distribution in anti-coincidence with fission, 
i.e. for particles emitted from nuclei which subsequently end up as evaporation residues. 
The initial compound nucleus is 7o°Yb at 293 MeV excitation energy and an angular 
momentum of 52 Ti. The curves for different types of emitted particles are displaced 
against each other because charged particles gain energy in the Coulomb field of the 
residual nucleus. Furthermore, the thresholds for emission of n, p, and a are in general 
different. 
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In addition, for a-particles, the spectral distribution for particles emitted in coin- 
cidence with fission is slightly shifted towards smaller energies as compared to the distri- 
bution obtained when measured in anti-coincidence with fission. This is probably due to 
the fact that charged particles are preferentially emitted from the pole tips around the 
long half axis. The larger deformation then implies a smaller gain of kinetic energy from 
the repulsive Coulomb field. 

In Figs. 20 and 21, the normalized yield for neutron emission is shown as a function 
of the nuclear deformation in the fission and evaporation channels. The reader should 
notice that for nuclei which undergo fission, the emission of neutrons takes place, on 
the average, at a larger deformation than for the nuclei which end up as evaporation 
residues. This is due to the fact that the distribution of the fissioning nuclei moves 
towards the saddle point in the deformation landscape. This effect should give rise to 
an experimentally observed anisotropy in the angular distribution of prefission particles 
different from the one observed in the angular distribution of particles emitted by the 
evaporation residua. The distribution for the case when the emission of protons and 
alpha particles is inhibited is shown in Fig. 20 by the short-dashed line. It is seen that 
this distribution is very close to the one obtained when the emission of all three kinds of 
particles is allowed. The initial angular momentum is different for the two Figs. 20 and 21. 
It is seen that the maximum of the distribution of tne neutrons emitted by the evaporation 
residua is shifted with increasing angular momentum towards larger deformations. The 
difference in the average deformation of the both distributions should be directly reflected 
in the difference in the anisotropy of the angular distributions of prefission neutrons and 
those emitted by the residua. 

3.5 Dependence of the multiplicities on the initial angular mo- 
mentum of the compound nucleus 

The multiplicity of prefission neutrons, protons and alpha particles is plotted in Figs. 
22-23 as a function of the initial angular momentum L. The initial compound nucleus 
is 160 Yb with the initial excitation energy E*=251 MeV in Fig. 22 and 293 MeV in Fig. 
23. The corresponding fission rates are plotted in Fig. 24. It is seen in Figs. 22 - 23 
that the neutron multiplicity decreases significantly with growing L while that for protons 
and alphas is much less affected. This is due to the fact that for large angular momenta 
the fission barriers Ub become small (see Fig. 25) and it takes a shorter time to reach 
the scission configuration. Consequently less neutrons are emitted on the average when 
the barrier is small. Alpha particles and protons are mostly emitted in the initial stage 
when the excitation energy of the nucleus is large, so that their multiplicities depend more 
weakly on the initial L. The effective fission rate (including particle emission) changes by 
2 orders of magnitude (see Fig. 24) when one increases the value of the initial angular 
momentum by 12% units up to its maximal value which corresponds to Ub=0- Fig. 23 
shows that this also holds true for the higher excitation energy of E* =293 MeV. 

This result indicates that a more precise knowledge of the initial spin distribution 
in compound nuclei is necessary in order to compare to the experiment. One has to stress 
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here that large L values contribute most to the experimental value of the multiplicities 
as already seen in Fig. 24. It therefore is certainly not correct to describe the experiment 
with a theoretical calculation using only the average value of the initial spin. As the 
fission barrier heights decrease strongly with increasing angular momentum as seen in 
Fig. 25, the higher L components which would lead to fission will have to be mocked up 
by an overestimation of the role of the fluctuating forces. 

Assuming an equal population of the initial angular momenta, which is, as just 
mentioned, not the best approximation, we have estimated the average of the n , p and 
a-multiplicities. The resulting values are compared in Table 1 with the experimental data 
taken from Ref. |J. All parameters of the model are those given in Sec. 2.1 and in Ref. 
|2"I|| . We only have to assume the preformation factor f a = 0.2 (see discussion in Sec. 
2.4) in order to reproduce the experimental number of alpha particles. This goes into the 
right direction since our calculations show that a particle emission is strongly enhanced 
by rotation and deformation effects. 

3.6 Influence of the friction forces on multiplicities of prescis- 
sion particles 

The role of the friction forces in fission of high excited compound nuclei is widely known. 
It was shown already in Ref. [f|] that the nuclear viscosity influence significantly the 
prescission neutron multiplicities. This effect was discussed later in almost all papers 
dealing with the problem. It was even assumed that the neutron multiplicities could be 
an indirect measure of the nuclear friction (see e.g. Ref. |§). 

The dependence of n , p and a-multiplicities on the magnitude of the friction co- 
efficient connected with the elongation of the nucleus is plotted in Fig. 26 for 160 Yb at 
E*=293 MeV and L=50h. The strength of the friction coefficient is varied in the interval 



(0.02, 2) in units of its value given by the wall formula fll4]| . The neutron multiplicity 
for small values of the friction coefficient is about 30% smaller than obtained with the 
wall friction. As the transient time increases as a function of the friction, the increase of 
the multiplicity of neutrons emitted prior to fission is easly understood. Again the effect 
upon the emission of p and a-particles is small, because these particles are emitted in 
a relatively short time interval after the formation of the compound nucleus which is in 
our case small compared to the transient time. One has to notice that the results for 
very weak friction depend significantly both on the width of the initial distribution of the 
momentum conjugate with the fission coordinate (eq. (|15|) ) and on the initial angular 
momentum. For larger values of the friction strength, close to jwaii, the influence of the 
initial distribution of the momentum p on the neutron multiplicity is weak (see also the 
discussion in Ref. [|l|]). 
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4 Summary and discussion 



Our results demonstrate the importance of nuclear deformation on the evaporation of 
light particles from strongly excited nuclei. This dependence on the deformation plays 
an important role also for the competition between fission and particle emission and 
might modify the limits which were determined for the nuclear friction force Q from the 
experimental data [6] on evaporation and fission. The dependence of the evaporation 
width on the rotational angular momentum was found to be negligible for the emission of 
n and p and rather small for the emission of a-particles. Let us, however, remind that the 
fission probability depends very sensitively on the angular momentum, as we have already 
emphasized in Ref. ]T|. Due to this strong dependence of the fission probability on the 
initial angular momentum, it becomes very important to get precise information on the 
angular momentum distribution of the initial ensemble of compound nuclei. The outcome 
of the competition between light particle emission and fission may depend strongly on the 
initial angular momenta. Consequently, we have to try to obtain theoretical information 
on the angular momenta of the initial nuclei by treating the fission process dynamically 

MM- 

Agreement between calculated and measured emission probabilities for a-particles 
can only be obtained, if an empirical preformation probability for ct-particles of about 
0.2 is assumed. One of the most important further improvements of the theory will be 
to evaluate this preformation factor within the temperature-dependent Thomas-Fermi 
approximation which underlies our theory [J27| . 

Another aspect which we intend to investigate is the angular dependence of emitted 
light particles when the angular momenta of the initial compound nuclei are aligned or 
the nuclei are polarized. In these cases it is conceivable that the angular dependence 
of the emitted particles shows a more pronounced dependence on the rotational angular 
momentum than the integrated yields we calculate in the present paper. One might also 
hope that the angular distribution of emitted particles depends on the deformation of the 
source nuclei sufficiently sensitively so as to determine the deformation from such mea- 
surements. Of course, experimental data on the angular distribution of emitted neutrons, 
protons, and a-particles from aligned rotating deformed nuclei would be of great interest 
for these studies. 
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Appendix: Evaluation of the transmission coefficients w v 

The explicit calculation of the transmission rates given by eq. (§) requires the 
knowledge of the average transmission coefficient w u (e,£; x) for the emission of light par- 
ticles from the nucleus into the continuum. In principle, a rigorous treatment would 
require the knowledge of w u (e,£;x) which appears in eq. (|) and should be calculated 
by averaging the transmission coefficients corresponding to the transmission of a particle 
emitted in a given direction from each point of the surface of the deformed, excited and 
rotating nucleus. This can be done in principle but is hardly possible in practice where 
these calculations must be carried out at each time step and over a very large ensemble 
of trajectories. We introduce the following simplified procedure which follows four steps: 

• We first calculate the transmission coefficients at three selected points (1,2,3) at 
which the main body-fixed axes (x,y,z) respectively cross the nuclear surface (see 
Fig. Al). The coordinate system is chosen such that the rotation axis coincides 
with the x-axis and the z-axis is the symmetry axis. 

For all possible values of the angular momentum of the emitted particle 

£ = \ fi x pi\ = r, L p i[{ , (Al) 

where pj., is the component of the vector p parallel to the nuclear surface, and for 
each of the points (i=l,2,3) w is calculated using the Hill- Wheeler WKB expression 

w v (e,l,£ x ;i) = (l + exp[- 2 ^ ~ e) ]) _1 . (A2) 

Here e is the single-particle energy and the barrier height Eb corresponds to the 
maximum of the potential Vtot i n the direction r» ± perpendicular to the surface. The 
potential in which the particle moves is given by: 

Vtot = V nuc l + V cent + VqouI ~ ^£ x (A3) 

i.e. by the sum of nuclear centrifugal and Coulomb (in the case of charged 
particles) single-particle potentials. The quantity hui appearing in (A2) is given by 

V m u 

• Once w u (e, £, £ x ; i) is fixed at the three points (i=l,2,3) one averages in a second step 
over all possible values of the x-component of £ which is the component along the 
rotation axis. Such a procedure seems to be reasonable once the angular momentum 
of the emitted particles is not measured experimentally. 

Through this average procedure one obtains 

e 

E w u (e,£,£ x ;i) 

w v {e,£-i) = 2£+l , z = 1,2,3. (A5) 
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• The transmission coefficients w„(e, l\ i) being determined at the three selected points 
(i= 1 ,2,3) one carries out in a third step an interpolation which allows to calculate 
the transmission coefficients to any point {9, ip) of the nuclear surface 

w u (e, l\ 9, p) = sin 2 9 ■ (wicos 2 p + w 2 sin 2 ip) + w 3 cos 2 9 . (A6) 

This interpolation of quadrupole type ensures that the transmission coefficient in 
the directions along the main axes are the same as those for emission in the opposite 
direction. 

• Finally one averages the coefficients over the whole nuclear surface: 

/ w u (e, £] 9, p) da 
Me, l\ X) = ~ • (At) 

s 

The procedure described above incorporates in an approximate way the effects of nuclear 
deformation and rotation on the transmission coefficients. 
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Table caption 

1. Results of prescission neutron, proton and alpha multiplicity model calculations for 
160 Yb at excitation energies 251 MeV and 293 MeV compared with the experimental 
data taken from Ref. || .The theoretical results are averaged over all initial angular 
momenta with a weight proportional to the corresponding fission rates. 



Figures captions 

1. Emission widths for neutrons (n), protons (p) and alpha particles (a) emitted from 
the hot, rotating compound nucleus 160 Yb (E*=50 MeV, L=40fr) as a function of 
the elongation. 

2. The same as in Fig. 1 but for E* = 150 MeV. 

3. The same as in Fig. 1 but for E*=250 MeV. 

4. Emission widths for neutrons emitted from 160 Yb at E*=150 MeV and L=0, 20, 40 
and 60h as a function of the elongation. 

5. The same as in Fig. 4 but for protons. 

6. The same as in Fig. 4 but for alpha particles. 

7. Emission widths for neutrons emitted from different isotopes of Gd and Yb at 
E* = 150 MeV and L=A0h as a function of the elongation. 

8. The same as in Fig. 7 but for protons. 

9. The same as in Fig. 7 but for alpha particles. 

10. Probabilities to emit a neutron with energy smaller than e n from the fissioning 
nucleus 160 Yb for three deformations corresponding to the spherical shape, the top 
of the fission barrier and the scission point. 

11. The same as in Fig. 10 but for protons. 

12. The same as in Fig. 10 but for alpha particles. 

13. Probability to emit a neutron with energy smaller than e n for three excitation en- 
ergies (E*=50, 150 and 250 MeV) of the fissioning nucleus 160 Yb. 

14. The same as in Fig. 13 but for protons. 

15. The same as in Fig. 13 but for alpha particles. 
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16. Number of trajectories out of sample of 10.000 leading to fission as a function of 
time with particle emission turned off. The calculation is done for three different 
temperatures (T=3, 4 and 5 MeV) of the fissioning nucleus 160 Yb. The heights of 
the corresponding fission barriers (Ub) are indicated in the figure. 

17. Multiplicity of the prefission neutrons (n), protons (p) and alpha particles (a) as a 
function of time. 

18. Ratio between the number of trajectories leading to fission and the total number of 
trajectories as a function of time for 160 Yb. The emission of the light particles is 
taken into account. 

19. Energy spectra of neutrons (n), protons (p) and alpha particles emitted in coin- 
cidence with the fission events (solid lines). The dashed lines correspond to the 
spectra of the particles emitted from nuclei which end up as evaporation residua. 

20. Yield of neutron emission P n as a function of the deformation of the fissioning nucleus 
160 Yb with initial excitation energy E*=293 MeV and initial angular momentum 
L=52h. The solid line corresponds to the trajectories which lead to fission while the 
dashed line to those leading to evaporation residua. The dotted line describes the 
distribution when the emission of protons and alphas is turned off. 

21. The same as in Fig. 20 but for L=A5h. 

22. Multiplicity of prefission neutrons (n), protons (p) and alpha particles (a) as a 
function of angular momentum of the compound nucleus 160 Yb with the initial 
excitation energy E*=251 MeV. 

23. The same as in Fig. 22 but for the excitation energy E*=293 MeV. 

24. Ratio between the number of trajectories leading to fission and the total number of 
trajectories as a function of the initial angular momentum for two different excitation 
energies E*=251 MeV and 293 MeV. 

25. Fission barrier heights as a function of the initial angular momentum for two differ- 
ent excitation energies E*=251 MeV and 293 MeV. 

26. Multiplicities of prescission neutrons (n), protons (p) and alpha particles (a) as a 
function of the strength of the friction force. 
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Table 1 





E*=2hl MeV 


£*=293 MeV 


V 


model 


exp. 


model 


exp. 


n 


5.98 


6.10± 1.5 


7.80 


8.50± 1.6 


P 


0.94 


0.51± 0.07 


1.19 


0.70± 0.08 


a 


0.58 


0.48± 0.07 


0.66 


0.75± 0.08 
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Figure 2 
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Figure 4 
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Figure 7 
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Figure 17 
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